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ABSTRACT 



^S| ■ We use large volume, high resolution, N-body simulations of 3 different ACDM 

> ■ 

models, with different clustering strengths, to generate dark matter halo merg- 



ing histories. Over the reliable range of halo masses, roughly galaxy groups 



Q^ \ to rich clusters of galaxies, we quantify the number density of major mergers 

^^ . for two different time intervals and compare with analytic predictions based 

O 

I ' Key words: cosmology:theory - methodsmumerical -dark matter: merging 

o- ,. . 

$-H ' histories 

•^ : 1 INTRODUCTION 

c^ ■ It is now widely accepted that the observed large-scale structure in the universe formed hier- 

archically, through the process of gravitational amplification of small initial perturbations in 
a universe with predominantly cold dark matter (CDM). Within this "bottom- up" paradigm, 
larger structures form by the merging of smaller units, resulting in a dynamic and continu- 
ally evolving matter distribution. Multiple lines of evidence support this scenario, suggesting 
that clusters and groups are still forming in the present universe, and that the idealization 
of a relaxed, virialized structure is somewhat unrealistic. Thus the merger history is of fun- 
damental importance in determining not simply the details of structure formation but also 
in many cases the global properties of clustered objects under consideration. 

Among merger events there are two extreme limits: major mergers and accretion. In a 
major merger, the masses of the progenitor halos are comparable, and consequently their 
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interaction results in a halo which is dynamically disrupted for some time after the merger. 
In contrast accretion is the merger of a small halo with a much larger halo, which is generally 
not severely disrupted by the process. We shall focus throughout on major mergers of objects 
"in the field" (which are expected to differ from major mergers of sub-haloes within a larger 
halo, e.g. galaxies in a galaxy cluster, Cavaliere & Menci ( |I997| )). 

The dynamical disruption caused by a major merger is expected to have many observa- 
tional consequences (for an extensive list see Roettiger et al. ( |1996|) ). For mergers of groups 
and clusters, these include among other things: inducing scatter in IGM temperature ( |Math- 



iesen fc Evrard 2000|) , increasing and then decreasing star formation rates (e.g. Fujita et 



al (|1999|) , for galaxy size halo consequences see Barnes and Hernquist ( |1992|) , and refer- 
ences therein), disrupting cooling flows (e.g. Allen et al. (|1999| )) and producing non-thermal 
radiation ( |Blasi 200q ). 

For all of these reasons there has been intense study of mergers analytically, numerically 
and observationally. On the theoretical side these range from in depth studies of specific 
mergers ( [Huss, Jain fc Steinmetz 1999| ) to properties of mergers as an ensemble, their rates, 



consequences for formation times, etc., e.g. in Bond et al. ( |1991| ), Lacey & Cole ( |1993| ; 



T^MD , Kitayama & Suto (|TO|; |OT§), Tormen (lOTj), Somerville et al. ( |2DUq ), Percival & 



Miller ( |1999| ). This paper is an instance of the latter, that is, a study of the bulk statistical 
properties of major mergers in the field. Specifically the question we address is: What is the 
number density of objects, at a given time and mass, that have undergone a major merger 
within the last 0.5 or 2.5 Gyr? Part of our aim is to address this question with several 
different methods, to further understand our answers and to investigate the accuracy of 
some analytic estimates based on "extended" Press-Schechter (( |1974| ); hereafter PS) theory 
(see below). 

The outline of the paper is as follows. We discuss the simulations in §^ and the numerical 
merger counts in §^. A review of some of the relevant analytic work is given in §^ and these 
estimates are compared with our numerical results in §^. Some technical details are relegated 
to an Appendix. 

2 SIMULATIONS 

We have chosen to focus on high mass objects, roughly galaxy groups to cluster scales, 
where we believe that the evolution is dominated by gravity and thus relatively inexpensive 
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to compute numerically and under reasonable theoretical control. Our choice of time-scales 
is motivated by the fact that 0.5 Gyr is close to that relevant for star formation caused 
by mergers (e.g. as noted in Fujita et al ( |1999| ) and for galaxies argued by Percival & 
Miller ( pOOOp using luminosity curves in Bruzual & Chariot ( |1993| )) while the larger time in- 
terval, 2.5 Gyr, is closer to those expected for cluster relaxation processes (e.g. see Mathiesen 
& Evrard (gOOOD ). 



We ran three simulations of a ACDM model with a high resolution N-body code (Pagla 



19991) . The cosmological models all assumed ilm = 0.3 and ^a = 0.7 but differed in z = 
clustering strengths. We generated one realization for each of ag = 0.8, 1.0 and 1.2, 
where af is the variance of the matter fluctuations in top-hat spheres of radius 8h~^Mpc. 
Each simulation employed 128^ dark matter particles in a box of side 256 h~^Mpc, i.e. a 
particle mass of 6.7x IO^^/i^^Mq, with a Plummer force law with softening length 200/;,^^kpc. 
The initial conditions were generated by displacing particles from a regular grid, using 
the Zel'dovich approximation, with the initial redshift was chosen so that the maximum 
displacement was (slightly) smaller than the mean interparticle spacing due to the large 
volume of the simulation. Time steps were chosen to be a small fraction of the shortest 
dynamical time of any particle in the simulation. The initial redshift and number of steps 
ranged from Zi^^t = 9 and 445 time steps for erg = 0.8, to zi^^ = 14 and 635 time steps for 
(Ts = 1.2. Further details (transfer function used, etc.) are given in the Appendix. 

The large volume of the simulations provides good statistics on the high-mass end of 
the mass function, where we have focused our attention. We have checked that the resulting 
mass functions scale as expected over the range of masses and redshifts which shall be of 
interest in this work. In our merger counts we do notice some possible effects at early times 
for low mass halos in the erg = 0.8 simulation, and these would be consistent with residuals 
from our grid initial conditions and limited numerical resolution. We shall comment further 
on this in §^. 

2.1 Group catalogues 

The full particle distribution was dumped every 0.5 Gyr starting from z = 2.16. For each 
output we generated a group catalogue using one of two group finders. We used HOP ( [Eisen- 



stein fc Hut 1998| ) for all three runs. For comparison, as well as to make explicit contact with 
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other work, we additionally used the friends-of- friends algorithm (FOFQ, Davis et al. ( |1985| )) 
for the o"8 = 1 run. From these catalogues we were able to construct merger trees (see below) 
back to z = 1.86 for 0.5 Gyr intervals and back to redshift z = 1.11 for 2.5 Gyr intervals. 

Since it may be unfamiliar, we briefly review the operation of HOP. For technical details 
the reader is referred to Eisenstein & Hut ([19981). HOP finds groups by first assigning each 



particle a density, and then "hopping" to neighbors with a higher density. Each particle 
belongs to the same class as its densest neighbor, and in this way each particle is assigned 
to a local density peak. To correct for the possibility of local density maxima causing groups 
to fragment, groups are merged if the bridge between them exceeds some chosen density 
threshold. There are 6 parameters one must choose for HOP: we used the default values 
for A^merge and A^hop (number of neighbors to look at when searching for a boundary and 
the densest neighbor) but took a lower value of the number of particles averaged over when 
calculating the density (ri^ = 8 rather than 64). The results were similar for n^ = 8 and 
n^ = 16, with the higher Ud excluding small groups. We excluded all groups of fewer than 8 
particles. The remaining three parameters are all density thresholds. Eisenstein & Hut ( |1998| ) 
claim the method works best if the thresholds are in the ratio 1 : 2.5 : 3, so we took 
(^outcr = 50, (^saddle = 125 and 5peak = 150. Thcsc thresholds correspond to: the required 
density for a particle to be in a group, the minimum boundary density between two groups 
for them to be merged, and the minimum central density for a group to be independently 
viable. We experimented with these thresholds and found negligible differences in the final 
catalogues. At redshift zero there were between 18,000 and 19,000 groups (see Fig. |ID. 

The Friends-of- Friends algorithm is much better known. It has only a single free parame- 
ter: the linking length b, usually given in units of the mean interparticle spacing. EOF defines 
groups of particles which are each separated by less than b from at least one other member 
of the group. Roughly speaking a EOF group consists of all particles within an iso-density 
contour of b~^ of the background matter density. We follow convention and choose b = 0.2 
(corresponding to a local overdensity 125). For the as = 1 case studied here, some notable 
differences for merger counts and number counts were found, shown in Fig. |^. For halos with 
more than 50 EOF particles, EOF halos had on average about 90 per cent of the particles 
of the corresponding HOP halos. EOF did find a larger total number halos, but these extra 
halos were at very low mass. 



The implementation used here is from pttp://www- hpcc.astro.washington.edu/tools/FOF/ 
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Both of these group finders do not identify substructure within large halos. This is not a 
particular problem for our purposes since we are interested primarily in mergers which take 
place in the field. Mergers within clustered environments are expected to behave differently 
( Pavaliere fc Menci 1997| ) and have been studied by other groups, e.g. Gottloeber, Klypin 



& Kravtsov ( pPUq ) 



2.2 Merger trees 

The definition of a major merger is somewhat arbitrary as there is some dependence on what 
physical consequence is of interest. In this paper we define major mergers to be mergers in 
which the portions of the two largest halos that become part of the final halo have a mass 
ratio between 1:1 and 1:5.[J] We compute the number of mergers within a given time interval, 
a very interesting quantity from an observational and theoretical standpoint. Note that this 
is an inclusive rate - the number of mergers is not given directly by a merger rate since 
in a given time interval, halos can merge and then accrete, or accrete, merge, accrete, or 
merely merge, etc. We count any halo which has a major merger in the given time interval, 
regardless of what else occurred to it during the interval. 

Having obtained the group catalogues for each output we constructed a merger tree 
for each simulation. For each time and for each (parent) group, we identified the (daughter) 
group membership of all the constituent particles at the next output time. The parent group 
was then considered a predecessor of the daughter group if more than 8 parent-group particles 
belonged to the daughter group at the next output time. For 0.5 Gyr merger trees, four of 
these "daughter groups" were considered for each predecessor group. This resulted in missing 
very few of the daughter groups - less than twenty total daughter groups went uncounted 
for the combined 22 outputs times three values of ag for HOP (there were ^ 10,000 groups 
for most time steps). For the 2.5 Gyr merger trees, a more conservative maximum of 15 
daughters was recorded for each predecessor group and no predecessor group had more than 
this number. 



T Note that this definition also allows an unbound collection of particles from a parent halo to count as a predecessor for the 
purposes of defining a merger. Quantifying how closely this criterion restricts one to a resulting dynamically disrupted halo 
would be very interesting. We thank the referee for emphasizing this issue. 
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3 NUMERICAL MERGER COUNTS 

The runs and merger tree construction described above resulted in halo evolution and merger 
counts with the following properties. We found, in agreement with Tormen ( |1998| ), that about 



20 per cent of the mass of the halos was lost as the halos evolve. In fact some groups, in 
particular many of those close to the minimum group size, disappeared completely between 
time steps. While this could be partly due to our finite force resolution, it also highlights 
the ambiguity inherent in halo identifications using these group finders. 

We anticipate that small mass predecessors may be undercounted, especially at earlier 
times. For this reason we only present results for halos with more than 10 times the minimum 
group size, or 5 x 10^^ h~^MQ. f] When predecessor sizes are compared for the purpose of 
identifying a major merger (i.e. a ratio 1 : 1 to 1 : 5), the number of particles from the 
predecessor group that go into the final group, not the number of total particles in the 
predecessor group are considered. 

The number density of halos and of merger counts (within 0.5 Gyrs) at redshift zero for 
the three values of ag (using HOP), are given in Fig. ^ The smooth curve going through 
the upper set of points (the number density) is the PS fit discussed in the next section. As 
we can see, the PS theory provides an adequate fit to the simulation results, as has been 
noted before by numerous authors ( [Efstathiou et al. 1988| ; [Efstathiou fc Rees 1988| ; [White, 



Efstathiou fc Frenk 1993t |Lacey fc Cole 1994| ; |Eke, Cole fc Frenk 199^ ; [Gross et al.T998| ; 



Governato et al. 1999[) . The lower curves correspond to two of the analytic fits, also discussed 

in the next section, as well as a quadratic fit to logM/{h~^MQ) and ln(l + z) for all times. 

Many groups have found dependence on group finder of their results, e.g. comparing 



SO(178) and FOF with different linking lengths (e.g. Lacey and Cole ([1991 ), Tormen ([T998[) ) 



Our FOF merger trees differed from those for halos identified with HOP, just as the total 
numbers of halos differed. The comparison between the HOP and FOF groups and merger 
trees can be seen in Fig. |^. 

One can consider the number of mergers within a given recent time interval, taken here 
to be 0.5 Gyr, as a function of redshift, as is shown in Fig. ^. One can see that the number 
of high mass halos decreases as one goes back to earlier times, and that simultaneously the 

■f This was because one could find in the 0.5 Gyr merger trees major mergers with a second predecessor as small as 1/10 the 
size of the final mass. For the 2.5 Gyr merger trees of order 20 of the second predecessors were smaller than this fraction for 
medium range masses, for each value of ag, summed over all times. 
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Figure 1. Number density of groups and those with recent (within 0.5 Gyr) major mergers at z = found by HOP, for three 
values of erg. The triangles are the simulations, the upper (smooth) curve is the PS prediction (Eq. ih with 5c = 1.48. The lower 
curves are quadratic fits as a function of log mass and ln(l + z). The dotted line is the fit to the simulation results, the dashed 
line corresponds to the fit to the "direct" calculation, Eq. (H), and the dot-dashed line corresponds to the "jump" estimate 
given in Eqs. (|lC|). One halo in the entire simulation volume corresponds to a density of 6 X 10~*(/i~^Mpc)^. 

fraction of total halos of any given mass which have undergone a major merger within a 
recent fixed time interval increases. 

A related quantity, the merger counts^ as a function of redshift at fixed mass, is plotted 
in Fig. 0. It appears that, for lower mass objects, the number of merger counts within 0.5 



The merger rate, rather than counts, as a function of redshift was studied by e.g. Percival & Miller (1999), Percival, Miller 



Sz Ballinger (1999), along with consequences such as star formation. 
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Figure 2. Number densities of halos and mergers for the 0.5 Gyr time interval, using two group finders, HOP and FOF, on 
the eg = 1 simulation. At left is redshift z = 0, at right is z = 0.53. The best fitting Sc in the PS formalism for the HOP halo 
number density is 5c = 1.48 and for the FOF halo number density is 5c = 1.64. 

Gyrs first increases with look-back time and then decreases. For higher mass objects the 
number merely decreases, which may be due to the lack of available high mass objects at 
earlier times. The trend is similar for the 2.5 Gyrs merger counts, however for ag = 1.2 and 
low mass {10^^'^h~^MQ) it appears that there is not yet any evidence for a high redshift 
decline in counts. 

In the next section some analytic estimates of number of major mergers within a given 
interval will be made using (extended) PS theory. 



4 ANALYTIC ESTIMATES 

There has been much theoretical work on major mergers and on mergers more generally. 
The most popular, and successful, semi-analytic formalism for predicting merger rates and 
associated quantities is PS theory and its extensions ( [Bond et al. 1991| ; [Bond fc Myers 19"96| ; 
Bower 1991| ; |Lacey fc GoleTgg3| ; ^.acey fc GoleT^U^ ; |Kauil'man fc White 1993| ; [Kitayama 



fc Suto I996a| ; [^omerville fc Kolatt 1999| ; [Ibrmen 1998| ; [^heth 1995| ; |Sheth 1998| ; |Sheth fc 
[Lemson 1999| ), although there have also been merger calculations done in other frameworks 
( Parlberg 1990| ; |Blain fc Longair 1993| ; [Cavaliere, Colafrancesco fc Menci 1992| ; Pavaliere fc 
[Menci 1997|) . Extended PS provides the necessary ingredients for calculating the quantity 
of interest here: the number of halos of a given mass which have had a major merger within 
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Figure 3. Number density of groups and major mergers within 0.5 Gyrs at z = 0.0, 0.53, 0.98, 1.42 for erg = 1. The curves are 
as in Fig. ll i.e., the upper (smooth) curves again are PS predictions for number densities with Sc = 1.48 and the lower curves 
are quadratic fitting functions (as a function of log mass and ln(l + z)) to merger counts in section §3: the dotted line is the 
fit to the simulation results, the dashed line corresnon ris t o the fit to the "direct" calculation, Eq. (^8h, and the dot-dashed line 
corresponds to the "jump" estimates given in Eqs. (|ia,hd). One halo in the entire simulation volume corresponds to a density 
of 6 X lO-^ih-'^Mpc)'-*. 



some given time frame at some given epoch, which is to be compared with the simulations 
described above. 

Press-Schechter and extended PS theory provide extremely useful characterizations of 
the distributions of mass and the growth of structure. They are computationally inexpen- 
sive and exhibit clearly the dependence on cosmological parameters. Their use has been 
widespread because, for many quantities (e.g. in Bond et al. ( p.991|) , Lacey & Cole ( [1993| ; 
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Figure 4. Simulated major merger counts as a function of redshift for 4 different masses. The number of recent mergers increase 
in time and then decrease, perhaps due to the fewer numbers of halos at earher times. The time period is 0.5 Gyrs for the left 
panel and 2.5 Gyrs for the right panel. 



, Kitayama and Suto ( |lDiJ(Ja| ; |1996b| ), Tormen ( [T^^ , Somerville et al. pnUD| )), they 



are in good agreement with numerical simulations. 

Given the simplistic assumptions going into the PS theory, it is intriguing that the 
agreement with numerical simulations is so good. Many of the assumptions are known to 
be incorrect in detail: spherical collapse QSheth, Mo fc Tormen 1999| ), the monotonic growth 
of halos (e.g. Tormen (|1998|) ), and the association of initial density peaks with final halos 
( lUarlberg 1990| ; |Katz, Quinn fc Gelb 1993| ; |Frenk et al. 1988|) . Hence, while PS is very 
useful and captures something essential about the process of structure formation, its failures 
( pross et al. 19981 |Sheth fc Tormen 1999| ; [Jenkins et al. 20"00| ) can tell us something very 
interesting as well. For example, extending the picture of spherical collapse to ellipsoidal 
collapse improves agreement with the mass function from numerical simulations ( |Monaco| 
[I997a| ; [Monaco 1997b| ; [Lee fc Shandarin 1998[ ; [Sheth, Mo fc 'Ibrmen 1999[ ). It is to be hoped 
that finding where and how extended PS theory fails to predict mergers will shed further 
light on this issue. 

We briefly review some background from the theory of PS and extended PS which we 
shall make use of below. More detailed definitions are in the Appendix. 
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4.1 Comoving number density of spherical collapsed systems 

The basic prediction of PS theory is the comoving number density of viriahzed halos with 
mass in the range (M, M + dM) at time t. The PS prediction is 



Nps{M,t)dM 



1 Po Scjt) 
/2^Ma3(M) 



da\M) 



dM 



exp 



m) 



2(t2(M) 



dM 



(1) 



Here po is the comoving, mean mass density of the universe, for matter density Qrn, ^c{t) 
is the threshold density contrast for spherical collapse at time t and a'^{M) is the vari- 
ance of the matter fluctuations smoothed over a region of radius R corresponding to mass 
M = 4nR^po/3. For the functional form of Sc{t) and a fit to cr'^iR) for the model under 
consideration see the Appendix. 

In order to compare the simulations to analytic PS estimates, Sc must be specified. The 
spherical top-hat collapse model (Peacock 1999| ; fiddle fc Lyth 2000| ) predicts that 6c — 1.69 
with a small cosmology dependence. We choose to regard 6c as a free parameter in the PS 
theory and adjust it to get the best agreement with the simulations for a given choice of 
group finder. The groups found by HOP resulted in a set of halo mass distributions fit within 
40 per cent (better than 20 per cent most of the time) by PS predictions for M > IO^^H^^Mq 
if 6c = 1.48 was taken. For FOF the best fit is obtained with 6c = 1.64, closer to the spherical 
top-hat prediction. 

4.2 Conditional probabilities 

The random walk model inherent in the PS theory can be extended to describe conditional 
probabilities: given a point in space that ends up in a halo of mass M2 > Mi at t2 > ti, the 
probability that it was in a halo of mass Mi at ti is: 



Pi(Mi,ti|M2,t2)dMi 



6r^ — 6. 



c2 



/2^(a2-ai)3/2 



daj 



dMi 



exp 



(5. 



cl 



Jc2 



2(af - al] 



dMi 



(2) 



where 6c2 = 6c{t2), 02 = o"(M2), etc. Given that a point starts in a halo of mass Mi at ti, 
the reverse probability that it ends up in a halo of mass M2 at t2 is 



P2{M2,t2\Mi,ti)dM2 



4.3 Merger rates 



1 Sc2i6cl - 6c2) 



/2n 



6. 



cl 



af 



c^K^i - ^2) 



3/2 



da^ 



dMo 



exp 



a|(5ci - af6\ 



2X2 
c2 



2a!al{aj - a^ 



dM2{3) 



Merger rates are calculated as a limiting case of the conditional probabilities above. The 
probability of a mass change over a small time is 
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Pi,2(M + AM,t^At|M,t) 



dPi,2(M^M + AM;t) 
Jt 



At 



(4) 



if At is small enough. (Note that Pi^2{M + AM, t|M, t) = 0.) One might think of a large 
mass change in a small time as coming from the addition of one halo, i.e. a merger (Lacey 
& Cole ( |1993| )). With this interpretation one gets the 

• Rate at which a point in a halo of mass Mi is incorporated into a halo of mass M2 at t: 



dPi{Mi->M2;t) 
dt 



dMi 



limAt_.o 

1 1 



Pi{Mi,t~At\M2,t) 
At 



dMi 



y/2^ {al- 



\ d5c{t)] 
dt 




da^ 
dMi 



dMi 



(5) 



• Rate at which a halo of mass M2 is formed from addition of mass to a halo of mass Mi 
at t, per point: 



dP2(Mi^M2;t) 



dt 



dM2 = limAt^o ^I^iMi^^MMi^ldM^ 



2lT 



\ "" 1 


3/2 


r dSc(t)] 

dt 




da^ 

dM2 


exp 


[a^{af-ai)\ 



Mtf 



1 

2^ 



dMo 



(6) 



These bulk probabilities can be combined to calculate the number of halos of a given 
final mass that have undergone a major merger (i.e. have predecessors with a given ratio) 
within a certain time interval, with the latter chosen to correspond to the relaxation process 
of interest. (We note that a modified version of the PS formalism that differentiates between 
(instantaneous) major mergers and accretion has been developed in ( [Raig, Gonzalez- Casado 
I&: Salvador-Solel998| ; [Salvador-Sole, Solanes fc Manriquel998|) as well.) 

Heuristically, one can count the major mergers which have occurred at a specific time 
for halos of a given final mass, whether or not they accreted mass before or after the merger, 
by taking the 

starting # of halos 

X merger prob. per unit time (7) 

X prob. to end up in final mass range . 
This is for a specific time and must be integrated to study a time interval. Plugging in the 
probabilities and integrating over time this becomes 



N, 



^^^^direct 



(M2, ti, t2)rfM2 = //f dt J d{Mi + AMi) / dMiNps{Mi,t) 



X 



dP2i 



(8) 

^^^(Mi -^ Ml + AMi;t)P2(M2,t2|Mi + AMi,t) . 

This will be referred to as the "direct" estimatcQI] In words, this is the probability that at 
some time t in a given time interval of size t2 — ti a starting halo mass of mass Mi increases 

'I After the bulk of the work in this note was completed, the authors became aware of a similar expression in jCavaliere, Mcnci 
fe Tozzi 19991). 
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its mass instantaneously by AMi, times the probability that the resultant Mi + AMi mass 
ends up, at t2, in a halo of mass of M2. Taking Mi to be the mass of the larger of the 
predecessors, the range of integration for AMi (the mass of the second largest predecessor) 
is determined by the chosen merger mass ratio. This allows for the possibility of a merger 
happening any time between ti and t2, with accretion before or afterwards, which is in some 
sense the key difference between this and the work of Carlberg (|1990|) as formulated by 
Lacey & Cole (|1993|) . Another differences is the range of integration for the masses taken. 



An equivalent way of writing the abov^ is 

//. p r/P ]\/f 

dtj d{Mi+AMi) J dMi-^iMi -^ Mi+AMi;t)Pi(Mi+AMi,t|M2,t2)A^P5(M2,t2)^(9) 

The limits on the integrals are as follows. The time integral is taken to correspond to 
the time interval of interest (0.5 or 2.5 Gyrs), and the range of AMi is determined by 
the mass ratios required for a major merger. So for our case, with ratios 1 : 1 to 1 : 5, 
Ml/5 < AMi <Min(M/ — Mi, Mi). For the dMi integral, the chosen mass ratio fixes the 
upper limit to be 5M2/6. The lower limit is not obvious, but it turns out that the dependence 
is rather weak. By considering the simulations one finds that Mi G [M2/3, 5M2/6] will cover 
most of the range of largest predecessor masses for 0.5 Gyrs. This was also used for the 
2.5 Gyr run. To use this formula more generally one would like to be able to choose this 
lower limit without relying upon simulations. Changing this lower limit for the 0.5 Gyr case 



between M2/2.5 and M2/3.5 produced indistinguishable results.^ To evaluate Eq. (|]), two 



of the three integrals can be done exactly, but the last integral over Mi seems to require 
numerical integration (see Appendix). A quadratic fit to the corresponding logA^mer(> M), 
for 6c = 1.48 as a function of log(M) and y = ln(l + z) is shown as the dashed line in Figs. |l| 
and ^ 

It should be noted that even within the PS description this expression is only approxi- 
mate. If a halo has two jumps before getting to the final mass and they are both a sizeable 
fraction if its mass, they may be counted twice. In principle one could bound this effect by 
calculating the probability of two jumps occurring using a generalization of Eq. (^). One will 
also double count the number of mergers where the two initial masses are identical, but this 

II The factor M2/M1 is due to the one to one correspondence (by assumption) between initial halos of mass Mi and final 
number of halos. If a halo of mass Afi ends up in a halo of mass A/2 then all the mass in the Mi halo ends up in the halo of 
mass A/2- So one can count the number of initial, and thus final, halos by taking the mass fraction and dividing by A/i to get 
the corresponding number of halos. 
** We thank the referee for suggesting we calculate this dependence. 
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14 Cohn, Bagla & White 

is expected to be a small number. Note that the distribution of mergers with a given mass 
change should be independent of the initial mass, using the random walk interpretation of 
the PS formulae ( [Percival, Miller fc Peacock 2000[ ). However, for a major merger the given 



mass change of interest depends on the mass of the initial halo, and in addition the total 
number of mergers depends on the number of initial halos (which depends on both time and 
mass), consequently some dependence on mass and era is expected. 

There are other approximate expressions which can be considered. Here we focus on one 
which has a particularly simple origin, in order to see how well it captures general trends 
and the quantity of interest. This is to use a jump to M instantaneously as an estimate of 
merger counts. Saying that the largest mass component of the two halo merger has at least 
half the final mass and allowing it to jump by a mass in the range (M/2, 5M/6) corresponds 
to taking 

iV__(M,t)^ /|^dM /^-(^V^^-^) fiVp,(M,t) (10) 

This will be referred to as the "jump" estimate. Requiring the largest predecessor to have 
mass only > M/2 avoids counting changes in both components of a major merger as distinct 
mergers (except in the case, probably rare, where both initial masses are exactly equal to half 
the final mass). However, there is some under- counting expected because, as noted earlier, 
the simulations showed that the largest predecessor could sometimes be as little as 1/3 the 
final mass. This expression only counts the masses which come exactly to a given mass at 
the given time, rather than any masses which might have merged and then later accreted to 
reach the given mass, or those which have reached this mass and then accreted out of the 
mass range. As the time dependence factors out in front, the only difference between 2.5 
Gyrs and 0.5 Gyrs is the multiplication by an overall factor of 5 to change time units. 

This turns out to be a simple quantity since for the range of interest the effect of the 
factor M/Mi in Eq. (|10|) gives just an overall pref actor of approximately 1.44: 

iV^er,_(M,t) ^1.44/|^'rfMi^^(^^^i^iVp,(M,t) 



M 
2'- . "- , qi] 



1.44 



1 



d&c{t) 



dt 



NpsiM.t). 



Another approximation which is sometimes used as the formation rate (e.g. Kitayama 
and Suto (|1996|))Q, 



TT There is no factor of M/Mi here because the number of initial halos is not in one to one correspondence with the number 
of final halos. For example, a huge number of infinitesimally sized halos can go to form one final halo. 
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iV^er,_(M, t; M/2) = jf/' dM^^^^^^^p^NpsiM, t) 



'^ .J^^^(MI'£y^^MjA) 



dScjt) 



dt 



Nps{M,t). 



(12) 



which turns out (as seen by numerical integration in the regime of interest) to also be close 
to the "jump" estimate when multiplied by 1.44: 

^^^^jump ' ^^^^ form ' \ ) 

The "jump" estimate is somewhat larger, by about 4 per cent, for masses ~ 5 x 
10^^/i^^Mq and ~9 per cent larger for masses ~ 1.5 x lO^^/i~^M0. In both of these, one 
is only looking at mergers going up to the final mass of interest, while in reality accretion 
should be taking halos both in and out of the range of interest. These quantities will turn 
out to work surprisingly well, perhaps because of some averaging out from the accretion 
effects. 

There has been related analytic work using extended PS to calculate "formation times." 
These differ from the merger counts under study here. A few differing formation time defi- 
nitions are used, e.g. when half the mass in the halo has the halo has assembled ( [Lacey fc 



Cole 199^ ; [Kitayama fc Suto 1996aD or when the assembly of mass in the halo ends ([Blain 



fc Longair 1993| ; pasaki 1994| ; [Percival fc Miller 1999| ). For example, a major merger by our 



definition could occur after half of the halo mass assembled, e.g. with a halo of half the final 
mass merging with a halo of one third the final mass. Analytic and numerical formation 
times have been compared, e.g. by Tormen (|1998|) . 

5 COMPARISONS 

The analytic formulae in section §3 for the number of mergers, the "direct" and "jump" 
calculations (Eqs. ^, |10|) in the previous section, were compared with the simulation results 



described in section §2 for mass starting at 5 x VS^^h~^ Mq and above. As mentioned earlier 
this minimum mass was chosen because a major merger could have a smaller predecessor 
of mass 1/10 the final mass, and so the minimum final mass was taken to be 10 times the 
minimum halo mass. The two analytic estimates differed from each other and consequently 
had different success in fitting the simulation results. The goodness of fit also depended on 
the time interval (0.5 or 2.5 Gyrs) considered and (very weakly) on the group finder used. 
Comparison with Eq. (^, the "direct" estimate, will be made first. The analytic expres- 
sions were integrated numerically and then fit quadratically (in log(M) and ln(l + 2;)) and 
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16 Cohn, Bagla & White 

these quadratic fits were compared with the simulations. A quadratic best fit to the simu- 
lations directly, for all masses and redshifts, could vary from the binned simulation values 
by up to 40 per cent in the 0.5 Gyr runs - thus exact agreement with any smooth curve is 
not expected to be better than this. For 0.5 Gyrs and using HOP, the quadratic fit to the 
"direct" formula was in broad agreement with that found in the simulations. The analytic 
calculation tended to over-predict the number of low mass (M < IO^^H'^Mq) halos which 
had major mergers, and this occurred more often at early times and lower erg, being within 
a factor of three in the worst case and within 40 per cent much of the time. This may be 
in part due to the known tendency of PS to overestimate the number of halos, and thus 
predecessors for major mergers, at low mass (Gross et al (|1998|) , Somerville et al. (|2000|) ). 
A competing effect is the "loss" of small halos in the simulations. Another factor may be 
the residual effects of the initial conditions, again strongest at early times and smallest Ug. 
For mass ~ lO^^M© the "direct" prediction was close, over-predicting for lower mass and 
under-predicting for higher mass. (For example, for erg = 0.8 and the 0.5 Gyr time frame, 
the predicted "direct" merger counts for final halos of mass M < 10^^ were too high by more 
than 40 per cent down to z ~ 1.) For the 2.5 Gyr runs, the agreement at low mass was a 
lot worse between the "direct" estimate and the simulations, with the number of mergers 
predicted for mass less that IO^'^H'^Mq over-predicted by 40 per cent for almost all times 
and values of ag. The fit improved with increasing erg. 

Surprisingly, the "jump" estimate (Eq. [l^) worked quite well. The "jump" estimate was 
unbiased for 0.5 Gyr but for 2.5 Gyr tended to overestimate the number of mergers. As 
noted, the "jump" estimate was slightly larger than 1.44 times the formation rate based 
estimate and closer to the simulations for 0.5 Gyrs. Note that the "jump" and formation 
rate based estimates have time dependence only in an overall factor out front, taking the 
time units to be the interval of interest. In particular, for the two intervals chosen here, 0.5 
Gyrs, and 2.5 Gyrs, the predictions only differ by a factor of 5. This scaling was not observed 
in the simulations, which differed by a smaller amount. Thus if the "jump" or formation rate 
based estimate works for one time, it will not work very well for the other. We found that 
it worked well for 0.5 Gyrs (i.e. within 40 per cent most of the time) but not for 2.5 Gyrs 
(consistently too high, centered at an overprediction of counts by ~40 per cent). 

The percentage deviation for the two time intervals and the "direct" and "jump" esti- 
mates, for (Tg = 1, is shown in Fig. |. Each line corresponds to a different time step, and 
horizontal lines mark 0, 20 and 40 per cent deviation from 0. Large spikes at high mass 
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Figure 5. "Direct" and "jump" pre(iicte(i merger counts vs. simulations: Eqs. (|q, |lO[) for 0.5 and 2.5 Gyrs. The liorizontal lines 
indicate 0, 20 and 40 per cent deviations between the simulation and the analytic estimate. The simulation counts have been 
smoothed over a range of .05 in log M before comparing with the analytic expressions. 

are often due to small number statistics (which are worse at early times). For the "jump" 
calculation, the trends for the fits (within 40 per cent for 0.5 Gyrs and centered at 40 per 
cent too high for 2.5 Gyrs) were the similar for the other values of a%. 

In Fig. 6, the RMS deviation of the calculated versus simulated counts and the average 
deviation are shown as a function of redshift for all three values of a% for both the "direct" 
and "jump" predictions. 

Comparing the HOP and FOF group finders for the (Tg = 1 case, some notable differences 
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Figure 6. RMS and average deviation as a function of redshift for erg = 0.8, 1.0, 1.2 for "direct" (top) and "jump" (bottom) 
fits, Eqs. (| |i|). 

for merger counts and number counts were noted, as mentioned earlier. However, the analytic 
estimates worked similarly for both group finders when the appropriate be was used. 

In conclusion, it seemed that both analytic estimates gave comparable numbers of merg- 
ers to what was found in the simulations. The "jump" estimate worked better than the 
"direct" estimate in matching the shape of the curve A^mcr(> ^), but over- predicts counts 
for all masses for longer times. It is of clear interest to pursue this, but doing so will require 
a larger mass and redshift range and larger simulations. 

It is also interesting to try to find some estimate of the merger fraction from the sini- 
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ulations alone, similar to "universal" fitting functions for numbers of lialos as functions of 
various parameter combinations which have been found to have a particularly simple form 
(e.g. Sheth-Tormen ( [T999D , Jenkins (pOOOD ). We found that, using z/ = 6c/cr{M,z)^ the 
combinations ^^^'I'-^g^^ ~ z/(l + zY'^^^^ for 0.5(2.5) Gyrs and ^%|^ ~ //(I + zY-^ seem 
to be promising possible scalings. As the relatively small number of mergers makes these 
quantities very noisy, bigger simulations and better statistics will allow further study of this 
question. 

6 CONCLUSIONS 

In this note we have concentrated on the question of how many halos of a given mass and 
time have had a major merger (ratio of two largest predecessors ranging from 1 : 1 to 
1 : 5) within the last 0.5 and 2.5 Gyrs, for the mass range 5 x 10^^h~^MQ ~ lO^^h'^MQ. 
Numerical simulations for 3 values of ag in a ACDM cosmology were compared with two 
analytical estimates. These counts can be related to observational processes with 0.5 Gyr 
and 2.5 Gyr relaxation times. For 0.5 Gyr, the "jump" analytic estimate, Eq. (|10|), was 
closest to the simulations, although the "direct" estimate, Eq. (^, was close for low enough 
redshift and high enough ag- Ratios of major merger counts were also considered. The main 
conclusion is that the PS estimates of merger counts, a very interesting quantity from the 
observational point of view, roughly agree with simulations. This suggests that future study 
of this approach may be fruitful. Many extensions and refinements are possible, such as 
using larger simulations, going back to higher redshifts, and varying parameters such as 
time intervals and mass ratios. The dynamical disruption caused by the major merger may 
be further characterized by including additional information about the predecessors and 
daughter halos. For instance, our criteria do not discriminate between the transfer of a 
bound subhalo or the stripping of outer (perhaps spatially very far apart and dynamically 
unrelated) particles of a parent halo to a daughter halo. The former would be expected to be 
more likely to cause more dynamical disruption.)^ In addition, the analytic calculations may 



be able to use the more accurate fitting functions for masses (Sheth and Tormen ( |1999 ) 



Jenkins et al ( pOOOp ), if these approaches can be generalized to extended PS conditional 
probabilities. 

■f-f In terms of the often used parameter M* , M = M* *-> u = 1. 

SS We thank the anonymous referee for suggesting we introduce this distinction. 
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APPENDIX A: PRESS-SCHECHTER AND EXTENDED 
PRESS-SCHECHTER 

Al Notation 

In numerical simulations it appears that many properties of the final density field are present 
in the initial conditions and are simply "sharpened" by the non-linear amplification of grav- 
ity. Press Schechter ( |1974| ) theory utilizes this in an essential way, and couples it with the 
theory of non-linear evolution of a uniform, spherical overdensity embedded in a homoge- 
neous universe. For this reason we need to be able to predict the scale- and time- dependence 
of density fluctuations in linear theory. 

We shall convert between mass and length scales using a top-hat filter in real space, 
i.e. the mass associated with a smoothing scale R is M = {Ati/2i)ttR^Pq where po is the 
comoving mean mass density of the Universe, 

po = 2.7755 X 10" fi„ {h-^MQ){h-^M^c)-^ . (Al) 

We assume that the initial fluctuations were Gaussian with zero mean. The variance of 
the fluctuations, smoothed over a region of mass M, is given by 

/■OO rjh 

a\M) = (x\M,z = 0)= — W\kR{M))A]{k) (A2) 

Jo k 

where W{kR) is the window function for top-hat filtering 
The power spectrum (today) is given by 






^sik) = ^dr^ = SUk)—T\k) (A4) 

^0 



which uses the primordial density spectrum dnik) ~ A;" ^, for n ~ 1. We take n = 1. For 
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the transfer function in both the calculations and the simulations we used (Efstathiou et 
al. (|99|)) 

T{k) = [l + {ak + {bkf^ + {ckfy] ~^'' (A5) 

where 

a = {QA/T)h-^Mpc, b = (3/r)/i~^Mpc, c = {l.7/T)h~^Mpc, v = 1.13 (A6) 

In practice cr^(M) was calculated up to an overall constant which is fixed by the choice of 
(jg. A good fit near R = 8/i"^Mpc is 



loglO I —2 
\ Or 



-1.4936 log i?8 - 0.46554(logi?8) - 0.0982478(logi?8)'' 



(A7) 



where Rg = {R/8h~^Mpc). This is a fit to numerical integration between 0.27 < log R < 1.49 
which is better than ~ 1 per cent. 

Collapsed halos are taken to be regions in the linear density field with density greater 
than some critical density contrast, 6c, when smoothed on a scale R. In practice we take 
account of the linear growth by holding the variance fixed and reducing the density threshold 

4(t(^)) = (l + ^)^^4. (AS) 



We found 6c{z = 0) = 1.48 was a good fit to the mass distribution of the halos found by 
HOP ( lEisenstein fc Hut 1998|) and 6c{z = 0) = 1.64 a good fit for those found with FOF 
( Pavis et al. 1985| ). The growth rate can be approximated as ( Parroll, Press fc Turner 1992| ; 
Viana fc Liddle 1999|) 



g{Q{z)) = -Q 



and 

Q{z) 



Qr 



70 140 140 



(1 + ^) 



1 -1 



'"l-fi„ + (l + 2)3n^ 

for a flat A universe. 



(A9) 



(AlO) 



A2 Formalism 

The Press-Schechter formalism has been extended to describe histories of halos (Bond et 
al. (|T99l , Bond & Myers (|1996D , Bower ([T99l| ), Lacey and Cole ( [19931) , dHH, Kitayama 
& Suto ( p96il) , Baugh et al. ([19981) , Somerville & Kolatt ([T999[) , Tormen ( [1998[) ). A useful 
description can also found in the textbooks by Peacock ( |1999|) and Liddle & Lyth ([2000| ). 
Press-Schechter and extended Press-Schechter theory have a heuristic random walk inter- 
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pretation based on consideration of trajectories. Each trajectory corresponds to the motion 
of a point in space under the addition of modes 5k (i.e. filtering the density field around 
each point in space with filters of steadily increasing k). The arguments are rigorous when 
this filtering corresponds to a top hat filter in fc-space, that is, adding 5k for larger and 
larger values of k. This produces a random walk for each point 5 = Y^o^k- One can relate 
k — i>R,M^ a and thus one talks about trajectories as a function of M or a. Those trajec- 
tories which have crossed some critical value oi 5 = 5c{t) at any given time are taken to be 
collapsed halos. In terms of the random walk picture, one can consider 5c as an absorbing 
barrier and the mass function as given by the distribution of trajectories in k as they are 
absorbed. The resulting number density is Eq. (|I]), 



Nps{M,t)dM 



1 Po 5c(t) 
/2^Ma3(M) 



da\M) 



dM 



exp 



m) 



2a2(M) 



dM 



(All) 



This is the number of halos, for the amount of mass in halos of this mass (i.e. the number 
of trajectories), one multiplies this quantity by M/po- 

This trajectory approach allows one to consider trajectories which have crossed the col- 
lapse threshold (absorbing barrier) at two different times, corresponding to being part of 
a halo of one mass at one time and another mass at another time, Pi(Mi, ti|M2, ^2) and 
P2(M2,t2|Mi,ti), inEqs. (0,0). 

These conditional probabilities can be multiplied together to get joint probabilities: 



Pi(Mi,ti|M2,t2)iVp5(M2,t2) X MzdM.dM, 



PO 



P2iM2,h\M,,t,)NpsiM,,t,) X Mi.dM,dM2 = PiM,,t,,M2,h)dM,dM2 
which obey,Q2] e.g., 

dM^ J dM2Pi{M^,ti\M2,t2)Pi{M2,t2\Ms,h) = P^{M,,ti\M^,ts)dM^ . 
( [Sethi 19941) . 



(A12) 



(AM) 



A3 Integration of merger rate 

The starting integral is 



™™ A crucial identity for this is 



d« e-"/-'-''"' = J-e-2^^ 



(A13) 
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(A15) 



//f dt J d{Mi + AMi) / dMiN{Mi,t) 
x^(Mi-.Mi + AMi;t)P2(M2,t2|Mi + AMi,t). 

Writing o"^(Mi + AMi) = af and integrating between times ti and ^2 corresponding to 

Sc{ti) = k one gets 



dMi daf f da. 



1 



1 



Ml dMi 7 27r (af - a2)3/2 (^^ - a|)i/2 
Consider the integral over af: 

do-? 1 I 



1 _e '("'-"2) 



(A16) 



/ 



2^ (<7^0?)3/2 (o?-<x^)V2 



1 



__(^1-«2)1 

e ^(^F^ 



(A17) 



[2:2 ^™ 



r[l 



'a;l/2(y_a,)3/2| 

where we have defined 
A = 5i — 82 ■, X = a} ■ 
We then have 

dx 



-AV{2x)l 



a. 



2 ! 



CTi 



CTo 



y^/y - X2 yy/y - xi 



For the second term we can further define z = 1/x, to get 



+ /;^M^e-^'f(yz-l)-3/2 



u3/2 



_41 
e 2!/ 



27rA(l-erf[^(^2-l/y)^/']) 



2e ^ ^ 
(^2-l/s/)i''2 

1/21 



^2 



^1)] 



(A18) 
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(A20) 
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